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Abstract 

Since Br adsorption on Au(lOO) displays an incommensurate ordered phase, a 
lattice-gas treatment of the adlayer configurations is not reliable. We therefore use 
density functional theory slab calculations to determine the parameters necessary 
for the construction of an off-lattice model. We compute and analyze the total energy 
and electron density as the lateral Br position and coverage are varied. This allows 
the calculation of the corrugation potential, the short-range lateral interactions, the 
dipole moment (long-range interactions), and the residence charge. From these pa- 
rameters, we construct an off-lattice model with no freely adjustable parameters. 
The simulation results compare remarkably well with experimental results. 
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1 Introduction 



In recent years, a variety of theoretical studies of electrochemical adlayer for- 
mation have emerged as a result of new in situ experimental techniques, which 
provide atomic-scale structural information about the adlayer. For a review 
of Monte Carlo studies, see Ref. [1]. Among these studies, hahde adsorption, 
in particular, has provided a wealth of fundamental insight into the structure 
and formation of electrochemical adlayers [2]. Bromine adsorption on Ag(lOO) 
stands out as one of the most extensively studied examples, both experimen- 
tally and theoretically [3,4,5,6,7,8,9,10,11,12,13,14] 

The Br/Ag(100) system displays a second order phase transition from a dis- 
ordered adlayer at low potentials to a c(2 x 2) adlayer at higher potentials. 
The most successful theoretical treatments of the Br/Ag(100) adlayer use a 
lattice- gas approximation [5,6,7], in which the adsorbates are assumed to oc- 
cupy an array or lattice of adsorption sites. However, when a system displays 
incommensurate ordered phases, the lattice-gas approximation is not a valid 
assumption, since the periodicity of an incommensurate phase and the peri- 
odicity of the substrate surface do not match. In these cases, an off-lattice 
model must be constructed, which consists of laterally interacting adsorbates 
in a two-dimensional corrugation potential. A full and formal description of an 
off-lattice model for Br/Ag(100) can be found in Ref. [11], and the lattice-gas 
approximation is seen to be a special case of the off-lattice model in which the 
amplitude of the corrugation potential goes to infinity. 

For a more complete understanding of simple adsorption phenomena, we should 
also consider systems which display incommensurate phases, and bromide ad- 
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sorption on Au(lOO) is a good candidate for studying commensurate to incom- 
mensurate phase transitions. In situ Surface X-ray Scattering (SXS) [15] and 
Scanning Tunneling Microscopy (STM) [16] measurements of Br adsorbed on 
Au(lOO) in solution, as well as ex situ Low Energy Electron Diffraction (LEED) 
measurements [17], confirm that Br adsorbs preferentially at the bridge site 
of the Au(lOO) substrate. Three phases are observed for Br on Au(lOO): a dis- 
ordered phase at low coverages, a commensurate ordered phase with coverage 
= 1/2, and a uni-axially compressed incommensurate phase with > 1/2. 
Two points of theoretical interest are immediately obvious. First, since the 
system displays the same binding site and ordered phases in both vacuum and 
in solution, the need to include water into the theoretical model is greatly 
reduced, as water does not seem to play an active role in the adlayer forma- 
tion. Second, the presence of an incommensurate phase excludes a lattice-gas 
treatment of the adlayer, in contrast to the Br/Ag(100) system. 

In previous theoretical studies, symmetry arguments and parameter fitting 
have been the primary methods used for deriving the lattice-gas model pa- 
rameters. These techniques have worked well for systems requiring only a 
single parameter, like Br/Ag(100) [5], but these techniques for deriving model 
parameters have met with varied success for more complex systems, like Cu 
and sulfate co-adsorption [18]. With the introduction of the corrugation poten- 
tial, necessary for an off-lattice model, the complexity of the parameter space 
is greatly increased, and therefore, in this paper, we use density functional 
theory (DFT) to estimate the off-lattice model parameters for Br/Au(100). 

Before discussing the details of the DFT slab calculations, it is important 
to point out several issues which may be of physical importance but which 
were not or could not be included in the present work. First, there is no 
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effect of an electrode potential in the DFT calculation methods used here, 
either as an external electric field or, as would be preferred, in the form of 
charging and discharging of the substrate. This efi^ect is not included in the 
calculations because it is difficult and time consuming, and there is not yet 
a generally accepted method for its inclusion [19]. DFT cluster calculations 
of a single halide adsorbate in the presence of an external electric field have, 
however, been performed by other authors [20], but since macroscopic scale 
experimental systems contain many adsorbates, we expect the coverage to also 
represent an important degree of freedom. 

The second issue of possible physical importance is the inclusion of liquid water 
into the supercell. This is currently possible using a hybrid DFT/Monte Carlo 
or DFT/molecular dynamics technique to perform proper thermal sampling, 
but this is extremely computationally intensive. Despite the lack of water in 
the supercell, the model parameters derived here seem quite reasonable, as 
justified by the results of the off-lattice Monte Carlo simulations, presented in 
Section 4. 



This paper is organized as follows. Section 2 discusses the details of the 
DFT calculations and the analysis methods used to interpret the DFT re- 
sults. Section 3 presents the off-lattice model parameters and discusses their 
reasonableness and physical significance. Section 4 presents the results of an 
off-lattice model for Br/Au(100) and compares these results to experiment. 
Finally, Sec. 5 presents a summary and conclusions. 
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2 Methods 



All DFT results presented here were calculated using the Vienna Ab Initio 
Simulation Package (VASP) [21,22,23] using a plane-wave basis set, a periodic 
slab geometry, the generalized gradient-corrected exchange-correlation func- 
tional (GGA, PW91) [24], Vanderbilt ultra-soft pseudopotentials [25], and a 
cut-off energy of 400 eV. All slabs contained three metal layers separated 
by six vacuum layers, each vacuum layer being the equivalent thickness of a 
metal layer. To examine possible coverage dependence, supercells of two dif- 
ferent sizes were examined, a 2(2.97) x 2(2.97) x 12.6 A cell (denoted as 2 x 2), 
containing 4 surface Au atoms, and a 3(2.97) x 3(2.97) x 12.6 A cell (denoted 
as 3 x 3), containing 9 surface Au atoms. The /c-point mesh was generated 
using the Monkhorst method with a 7 x 7 x 1 grid for the 2x2 cells and a 
5x5x1 grid for the 3x3 cells. Initial convergence tests indicated that these 
parameters are sufficient for the accuracy desired, as justified by the results 
in Sec. 4. 

The Au(lOO) surface normal was oriented in the +z direction. A typical DFT 
supercell is shown in Fig. 1. Unless otherwise indicated, all Au(lOO) slab atoms 
were located at their bulk positions with a face-centered-cubic lattice constant 
of 4.2 A, as determined by our own ah initio calculations. This yields a (100) 
surface lattice constant of a = 2.97 A, as compared with the measured value 
of 2.885 A [15], an error of 3%. 

In many of the most common applications of DFT methods, the positions of 
all or some of the atoms in the configuration are relaxed in such a way as to 
minimize the total energy. However, since we are interested in more quantities 
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than just the total energy and we wish to compare results for different cover- 
ages, a full relaxation of the initial atomic configuration is not desirable for two 
reasons. First, even though Au(lOO) is known to display a surface reconstruc- 
tion, the surface reconstruction is lifted by even moderate halide coverages 
[15], thus eliminating the need to reproduce such a transition in the calcula- 
tions by relaxation of the slab in the surface tangent directions. Second, early 
in this research, an examination of the charge density revealed very large sur- 
face dipole moments for buckled surface configurations, e.g. when the surface 
atoms were relaxed independently in the surface normal direction. These large 
surface dipole moments made analysis and understanding of the Br-Au(lOO) 
bond polarity extremely difficult and problematic. Therefore, when any relax- 
ations were performed, they were performed by keeping all slab atoms at their 
bulk positions, except for the top layer atoms (Br-containing surface layer) 
which were relaxed simultaneously in a plane, and the Br layer was relaxed 
in a similar way. We introduce the notations Ztop and zbt which indicate the 
z positions of the top slab layer and the Br adsorbate layer, respectively. As 
is well known, relaxation of the top metal layers in slab calculations gener- 
ally has a relatively small effect on the adsorption energy, the effect being of 
the order of 5 to 100 meV, which is in agreement with the energy differences 
observed in this work between the fully relaxed and unrelaxed substrate. 

In addition to obtaining total energies, DFT results also include the electronic 
charge density, Pe{x), which has the sign convention of being more positive 
where the density of electrons is higher. Figure 2 shows the electronic charge 
density, Pe{z) — J Pe{x) dxdy, for several binding site configurations. The 
location of the three metal layers and the Br adlayer are clearly visible, and 
outside of the surface region, few differences can be seen between the different 
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adsorbed states. 

2. 1 Charge Transfer Function 

Unfortunately, the electron density, pe{x) or Pe(^), when plotted in its unmod- 
ified form (see Fig. 2), is difficult to interpret, and little insight is gained from 
so simple a plot. Even the full charge density, including both electrons 

and ionic cores, gives little novel understanding, and only when the "back- 
ground" charge is eliminated, is the presence of the polar surface bond easily 
seen. We therefore introduce the charge transfer function, so called because it 
clearly shows where charge has been transferred as a result of the formation 
of the (polar) surface bond. The charge transfer function is defined as follows: 



where p(.T)ArBr-Au(ioo) is the full charge density of the adlaycr system with 
N adsorbed Br in the cell, Pj(x)Br is the full charge density of a single Br 
atom at the same position as that in the Br-Au bonded system, i indexes the 
N adsorbed Br, p(x)au(ioo) is the charge density of the Au(lOO) slab with all 
atoms at the same positions as in the Br-Au bonded system, and the subscripts 
"e" denote electron only densities, having the sign convention that positive 
indicates greater electron density. The charge transfer function is calculated 
using only charge densities for the same size cells. 

Since all atoms are located at the same positions for each of the charge den- 
sities in Eq. 1, the ionic core charges cancel, leading to the second, equivalent 

7 




(1) 




Pe(^)iVBr-Au(100) " Ej Pei(^)Br " Pe(^)Au(100) } , 



expression in Eq. 1. Also, since we here consider only charge neutral cells, the 
charge transfer function must integrate to zero. Figure 3 shows ^p{z) for the 
same configurations shown in Fig. 2, and the charge polarization at the surface 
is easily seen. 

2.2 Residence Charge 

The net charge on the adsorbate is one of the simplest results to calculate 
from the charge transfer function. Unfortunately, there is no unique scheme 
to assign charge to each atom; however, in the simplest method, charge above 
the polarization center is assigned to the adsorbate, and the charge transfer 
function is integrated over this region. We define the residence charge as 



where the limits of integration are chosen between the polarization center and 
the chargeless region of the vacuum. Since the charge transfer function must 
integrate to zero for charge neutral systems, the charge residing on the slab 
must be equal and opposite to q. 

2.3 Surface Dipole Moment 

Motivated by the success of the lateral dipolc-dipolc repulsion in describ- 
ing the ordering behavior of Br electrodeposited onto Ag(lOO) [7], we calculate 
the dipole moment by 
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— * 

where Z is the center of the dipole, and the integral extends over the entire 
range of the charge distribution. 

From Eq. 3, it would appear that the dipole moment is dependent on the 
seemingly arbitrary center of the dipole; however, we note that when the 

— * 

charge density goes to zero at the limits of integration, -p is independent of Z 
[26]. Since all charge distributions considered in this paper go to zero in the 
vacuum region, we choose the limits of integration in the z direction to be 
[—3.2, 9.4] A, and the integration in the surface tangent directions obviously 
extend over the full x and y ranges of the cell. 

Prom symmetry considerations, the dipole moment must be of the form p = pz, 
at least for the symmetric [x^y] positions (bridge, hollow, and top sites), and 
we therefore introduce the surface dipole moment and define it as follows: 

s = p. z = J zp{z) nbt-Au(wo)/N dz 
^ J z {NAp{z) + p(^)au(ioo) + Ej Pi{z)Br} /N dz 

(4) 

= / zAp{z) dz + J zp{z) An{ioo)/N dz + Y.iI zpi{z)Br/N dz 
= / zAp{z) dz , 

where the factor of is included to correct for multiple Br-Au(lOO) dipoles 
in the cell. Of course, the charge distribution due to a single isolated Br has 
no intrinsic dipole moment and likewise for a symmetric isolated slab. Since 
these terms cannot contribute to the dipole moment, we are left with the final, 
simplified form of Eq. 4, and the sign convention is such that the dipole points 
from positive charge to negative charge. 
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3 DFT Results and Discussion 



In this section we present results from analyzing both the total energy and 
the charge transfer function. We examine both the lateral dependence and the 
coverage dependence using a variety of configurations. For the lateral depen- 
dence studies, the lateral and vertical position of a single Br atom in a 3 x 3 
cell were varied, and for simplicity, the slab atoms were fixed at their bulk 
positions. From the lateral dependence of the total energy, we can construct a 
sinusoidal approximation of the surface binding energy, which we call the cor- 
rugation potential. From the charge distribution, we can calculate the lateral 
dependence of the surface dipole moment and the residence charge. 

To examine the coverage dependence of the total energy and the charge distri- 
bution, only adlayer configurations with four-fold symmetry were considered, 
and all Br atoms were placed at the preferred bridge binding sites, see Fig. 4. 
In addition to the four cells shown in Fig. 4, two additional 2x2 cells were ex- 
amined with four and eight adsorbed Br, corresponding to coverages of = 1 
and 6 = 2, respectively. Although these highest two coverages are expected 
to be unphysical, their calculation was necessary to determine the short range 
excluded volume repulsion, needed for an off-lattice Monte Carlo model. 

3.1 Corrugation Potential 

We define the corrugation potential as the total energy of the adsorbate sys- 
tem, minimized with respect to ZBr, as a function of x and y. Experiments 
in both the in situ and ex situ environments [15,17] confirm that Br prefer- 
entially binds at the bridge site, indicating that the corrugation potential is 
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lowest at this site; however, this knowledge alone is not sufficient to determine 
the symmetry or amplitude of the corrugation potential. DFT calculations 
must therefore be used to determine the corrugation potential. 

For each x,y location (see Table 1) z^r was varied systematically. Figure 5 
shows the total energy as a function of z^r for several important y locations. 
The total energy minimum for each {x, y) position, along with the z value of 
the minimum are reported in Table 1. For clarity, the total energy of the bridge 
site minimum has been subtracted from each, and the Br sits just above the 
third metal layer which is centered at ^top = 4.2 A. 

A table of discrete results gives little impression of the surface symmetry, and 
since the off-lattice Monte Carlo simulations (see Sec. 4) require the corru- 
gation potential to be defined at all (x, y) points, we construct a sinusoidal 
approximation to the DFT results of Table 1. Certainly, many different ap- 
proximations can be made; however, we favor the following simple form which 
reproduces the DFT results extremely well around the bridge site, the region 
of most importance. 



where Ai = 274 meV and A2 = 234 meV. 

For comparison with the original DFT results, values of the sinusoidal ap- 
proximation are also listed in Table 1, and the full sinusoidal approximation 
is shown in Fig. 6. The differences between the DFT calculations and the si- 
nusoidal approximation are expected to have little qualitative effect on the 




(5) 



+A2 [cos (^) + cos (^)] /2 
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equilibrium behavior of the system, since these differences appear primarily 
at the low-probability binding positions. We also comment that the symmetry 
and amplitude of U {x, y) are consistent with the observed ordered phases and 
phase transitions [15], as will be made clear in Sec. 4. 

3.2 Lateral Interactions 

Assuming only two-body interactions, the effective lateral interactions between 
adsorbed Br can be determined by analyzing the total energy for configurations 
representing a range of coverages. Two things must be noted. First, to simplify 
the analysis, Zx,o^ —4.2 A and zbv =6.4 A in all calculations. Second, since 
we wish to compare results at different coverages, it is imperative that each 
adlayer should have the same square symmetry, thus changing the coverage 
has only the effect of changing the Br-Br nearest- neighbor separation, qbt, as 
illustrated in Fig. 4. 

To separate contributions to the total energy from contributions other than 
lateral interactions, we assume that the total energy for each adlayer configu- 
ration can be approximated as 

Etot = ^slab + A^^Br + NB + ^ , (6) 

where £'siab is the total energy of the blank Au(lOO) slab, Ebt is the energy of a 
single Br atom, N is the number of Br in the supercell, B is the binding energy 
of a single adsorbate to the substrate, and $ is the total lateral interaction 
energy between adsorbates. The values of the terms £^siab and Ebi can be 
determined from calculations including only the Au(lOO) slab and only a single 
Br atom, respectively. 
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The results of this calculation are shown in Fig. 7 and include results for the 
four coverages shown in Fig. 4 plus two additional cells with coverages © = 1 
and 2, corresponding to every other bridge site occupied and every bridge site 
occupied, respectively. These later two systems were included to examine the 
short range repulsion, which is a crucial parameter needed to construct an 
off-lattice Monte Carlo model. A non-linear fit to the data shown in the figure 
indicates 

$/A^ + B = 7 X lO^aif - 2560 , (7) 

with units of meV, and thus, B — —2560 meV. Note that the fit is primarily 
determined by the first three points, where the repulsion is largest, and this 
indicates that the repulsion represents the effective short-range ionic core 
repulsion. 

3.3 Surface Dipole Moment 

The polar surface bond can be measured by integrating over x and y to obtain 
Ap{z), shown for various binding sites in Fig. 3. All sites show a clear polar- 
ization of charge in the surface region and an oscillatory charge in the first 
two slab layers; however, the relative magnitude of these two behaviors is quite 
different between the four sites and is quite sensitive to site coordination. The 
quarter and bridge sites are seen to be qualitatively similar, both displaying 
a pronounced charge polarization in the surface region with moderate oscilla- 
tions in the slab surface region. The top and hollow sites, however, are quite 
different from each other and from the quarter and bridge sites, the top site 
having the largest oscillations in the slab surface region and the hollow site 



13 



having the smallest oscillations. 

Figure 9 shows the surface dipole moments for all considered configurations 
in the lateral dependence study. There is very httle direct dependence on the 
lateral position, but there is a strong, almost linear, dependence on zbt- The 
heavy hne in Fig. 9 indicates a hnear fit to all s values and has a slope of 
1.22 e, where e is the elementary charge unit. This suggests that Br/Au(100) 
forms an ionic bond, regardless of the coordination with the substrate. A two- 
dimensional slice through Ap, Fig. 8, confirms the p-type character of the 
electron-accepting orbital on Br. 

Although there is no direct dependence of s on the {x, y) position, the total 
energy does have a strong dependence on (x, y) , and this leads to an indi- 
rect lateral dependence of the surface dipole moment, when only the most 
energetically favorable z-qt are considered. These values are summarized in 
Table 1, and from these values, we can construct an approximate sinusoidal 
form needed for off-lattice Monte Carlo simulations, 

s{x, y) = 0.08 {h{y)g{x) + g{y)h{x)} + 0.14 , (8) 

where h{x) — cos(27rx/a) -|- 2, g{x) — cos(27rx/a) -|- 1, and the units are eA. 

For comparison with the DFT results, the values ol s{x,y) at several lateral 
positions are also listed in Table 1, and the sinusoidal approximation is ob- 
viously very good at positions with low values of the corrugation potential, 
where the adsorbate is preferentially located. For the top and surrounding 
sites, the approximation given by Eq. (8) is quite poor; however, since the 
adsorbate spends very little time at positions where the corrugation potential 
is large, the top region is not expected to significantly influence the behavior 



14 



of the adlayer. 

The surface dipole moment, s, can also be calculated as a function of coverage, 
and the results are shown in Fig. 10 for all values of ^top and zbt- The shaded 
symbols in the figure indicate the configurations corresponding to the energy 
minimum for each coverage. The results are summarized in Table 2. These 
results are to be expected from a simple examination of Ap{z) for the different 
coverages. See Fig. 11. No significant differences are seen between the four 
curves, suggesting that there is no significant coverage dependence of either 
the dipole moment or the residence charge. This can also be confirmed by the 
hnear dependence of the Fermi energy, Ep, on coverage. See Fig. 12. 

3.4 Residence Charge 

The lateral and coverage dependence of the residence charge are summarized 
in Tables 1 and 2. From these data, we can conclude that q — —0.12 ± 0.01, 
independent of coverage, with only bridge-site binding considered. Assuming 
that the double layer is dominated only by the Helmholtz layer, the residence 
charge can be related to the electrosorption valency as 7 = z — q/e = —1 — q/e, 
where z is the charge state of the ionic adsorbate species in solution ^ . This 
leads to a rough estimate of the electrosorption valency as 7 = — 0.9±0.1, 
with no direct coverage dependence, and this estimate still allows for the 
possibility of a complete, or nearly complete, discharge of an adsorbed Br~ ion 
as suggested in Ref . [15] . However, since these calculations do not include the 

^ By assuming a parallel plate capacitor of area. A, total charge, qQA, a plate 
separation, d, and a dipole moment s = qd, Eq. 18.18 of Ref. [27] reduces to 7 = 
z - q/e. 
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effect of excess surface charging or a external field, the connection between the 
calculated residence charge and the electrosorption valency is still unknown. 

A classical charge/image charge interpretation of this system is not correct. 
This can be illustrated more completely by examining the behavior of a point 
charge approximation of /S.p{z), which can always be constructed, even when 
the system is not charge/image charge like. The point charge approximation 
consists of two point charges of q and —q, located at the centers of charge of 
the adsorbate and the slab, respectively, and separated by a distance d, such 
that s = qd. As shown in Fig. 10, s has a linear dependence on zbt with slope 
ds/dzsr — 1.22e. If the charge/image charge interpretation is correct, then 
q{zBr) would be a constant, and the slope would be ds/dzBr = 2g. However, 
as Fig. 13 shows, q{zBr) is not a constant, clearly indicating that the system 
is not charge/image charge like. 

4 OfF-lattice Simulations 

The off-lattice Monte Carlo simulation method used here is discussed in com- 
plete detail in Ref. [11], and there is no need to repeat such a detailed discus- 
sion here. The simulated surface was an L x L square with periodic boundary 
conditions, where L was the number of surface atoms in a surface lattice di- 
rection. The corrugation potential used in the off-lattice simulations is given 
in Eq. 5, and the pair- wise lateral interactions between adsorbed Br are the 
summation of the strong short-range 1/r^ repulsions from Eq. 7 and a much 
weaker long-range dipolar repulsion. 

There are several different possible treatments for the dipolar repulsion. One 
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can assume that the dipole moment is flat or constant across the surface, 
that is s{x,y) = s, or one can assume the functional form of Eq. 8. Also, 
one can assume a simple dipole repulsion, known as the point dipole 
approximation, as used in Refs. [7] and [11], or alternately, one could use the 
more complex physical dipole form, involving the electrostatic repulsion of 
four different point charges. In this paper, we take the simplest assumption 
for the dipole-dipole repulsion, namely that of a constant dipole moment, 
s{x,y) = Sbridge = 0.30 eA, with the simple 1/r^ point dipole form used in 
Refs. [7] and [11]. The lateral interactions are truncated for r > 5a, and for 
r < 5a, 

(f)(r < 5a) = + Csr"^ , (9) 

where r is the distance between interacting adsorbed Br, measured in A, s = 
0.3 eA, and Ci = 14387 is the conversion factor which gives in units of meV, 
and C2 = 300000 is the constant which gives in units of meV. 

The Q{fl) isotherm results of the this model are shown in Fig. 14, along with 
four insets showing the two-point correlation function for each of the four 
observed phases. The two-point correlation function, C{Arx, /^Vy), is simply 
the probability of finding an adsorbate at a relative position of (Ar^., Ar^) from 
any other adsorbate. Thus, there is always a sharp peak at C(0, 0), which is 
the self correlation point. To improve the visualization, the self-correlation 
point, C(0, 0), has been removed from the insets, and the grayscale of each 
inset has been scaled between the highest value of C{Arx, Avy) (white) and 
the lowest value (black). 

The correlation functions in Fig. 14 can be viewed as a kind of average snap- 
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shot of the adlayer and are related to surface X-ray scattering data [7] . Prom 
left to right, the correlation functions in Fig. 14 represent a disordered com- 
mensurate phase {fi = 600 meV), an ordered c{y/2x 2-\/2)-R45° commensurate 
phase with Q — 1/2 (jl — 1000 meV), and two apparently different uni- 
axially compressed incommensurate phases (/i = 1060 and 1230 mcV), which 
represent a distortion between the commensurate ordered phase and an in- 
commensurate hexagonal phase. 

The four ordered phases observed in the simulations are quite similar to those 
found experimentally [15], except that experiments have observed only one 
incommensurate ordered phase. This apparent difference may signify that our 
model needs some minor adjustments, such as using a different approximation 
for the dipole-dipole interactions, or it may be simply an artifact of the square 
simulation geometry. The simulated surface is square with periodic boundary 
conditions, but the incommensurate phases are nearly hexagonal. This differ- 
ence in geometry obviously creates a strain on the incommensurate adlayer 
phases, and this strain may cause slight differences in the observed incommen- 
surate phases. The strain should decrease with increasing L; however, since 
these simulations for L — 32 required about six weeks of computational time 
on a PC with a single GHz processor, the effect of larger system sizes has not 
been checked. As such, the results shown in Fig. 14 should not be considered 
to be the definitive Monte Carlo results of this off-lattice model. 

The reader should note that off-lattice simulations are extremely computa- 
tionally intensive [11], and commensurate to incommensurate phase transi- 
tions depend very sensitively on the specific values of the parameters, thus 
making a reproduction of experimental results very difficult. Our preliminary 

models for Br/Au(100), based on less accurate Monte Carlo relaxations and 
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sampling, smaller system sizes {L — S and L — 16), and less accurate DFT 
parameter estimations, usually showed only one incommensurate phase. How- 
ever, we prefer to show the more accurate Monte Carlo results of Fig. 14, 
even though these may not fit as well to the experiments as the less accurate 
preliminary results. 

Figure 15 shows the experimental cyclic voltamogram (CV) from Ref. [15] com- 
pared to the simulated quasi-equilibrium CV. The simulated quasi-equilibrium 

CV current is 

.(^^_^de(^=B-^_ (10) 



where ^4 is a constant related to the area of the electrochemical cell and the 
scanrate, fl is the electrochemical potential (See Ref. [11].), B is a constant 
related to the choice of electrochemical reference potential, 7 = —0.9 is the 
electrosorption valency as determined from the residence charge, and E is the 
experimental electrode potential. 

The simulation curve is fit to the experimental curve in Fig. 15 by first calcu- 
lating dO/djj, directly from the Monte Carlo simulation. Then, the parameters 
A and B were varied, until the simulation and experiment curves matched near 
the phase transition from the disordered Br adlayer to the commensurate or- 
dered adlayer. Since the adjustment of A and B represents simply a scaling of 
the current and shifting of the potential, A and B should not be considered 
as true physical parameters of the off-lattice model, and as such, the results 
in Fig. 15 represent a model with no freely adjustable parameters, since all 
model parameters were estimated from DFT calculations. 
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5 Summciry and Conclusions 



Using DFT slab calculations, we have studied both the lateral and coverage 
dependence of Br adsorption on Au(lOO). We have examined both the total 
energies and the electron densities, and their analysis has enabled us to obtain 
the corrugation potential, the lateral interactions at short range, the dipole 
moment (the interactions at long range) , and the residence charge (the charge 
remaining on the adsorbate). 

We have found that the raw electron densities are difficult to interpret, and we 
have therefore defined and analyzed the charge transfer function, from which 
we can calculate the surface dipole moment and the residence charge. Both the 
surface dipole moment and the residence charge were found to show a strong 
dependence on the lateral Br position, but neither showed any noticeable 
coverage dependence. 

It should be noted that the coverage dependence study does not include the 
effect of the electrode potential or external electric field, and thus, a full un- 
derstanding of the coverage dependence, or lack thereof, is not yet known. We 
can, however, conclude that the lack of coverage dependence in the electron 
densities is not particularly surprising when one considers the energy scales of 
the adsorption problem. Electron binding energies are typically of the order of 
eV, whereas the lateral interaction energies at half coverage are of the order of 
tens of meV. From these energy considerations, we can conclude that simply 
changing the coverage has httle effect on the electron density. 

Considering the number of parameters and the sensitivity of off-lattice simu- 
lations to the values of the parameters, we feel that the comparison between 



20 



simulation and experiment is quite remarkable, considering that a number of 
presumably important effects were left out, such as an external electric field 
and the presence of liquid water. We also feel that the similarity of the simula- 
tion and experiment justify the values of the DFT parameters, and despite the 
differences between the simulation and the experiment in the incommensurate 
phase, constructing an off-lattice model from DFT parameters seems to be a 
valid and predictive method. Whether this technique will prove valuable for 
other systems remains to be seen. 

Further studies should include the effects of both water and the electrode po- 
tential to give further insight into these effects. However, since Br/Au(100) 
displays the same ordered phases and binding sites in both vacuum and aque- 
ous environments [15,17] and since our simulation results in vacuum compare 
so well to experiment, the importance of including water in this DFT study 
may be less important than in systems which are known to be strongly sol- 
vated, such as with systems containing sulphate. 
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Table 1 



Summary of results for the Br lateral dependence study. The nearest Au surface 
atoms are located at {x/a,y/a) = (0,0), (1,0), (1,1), and (0,1). The DFT results 
are given along with the values of the analytic approximations. 



x/a 


y/a 


ZBr 


Name 


Energy [meV] 


U{x,y) 


s [eA] 


s{x,y) 








[A] 




from DFT 


[meV] 


from DFT 


[eA] 


from DFT 








6.7 


top 


257 


508 


0.68 


1.10 


-0.19 





0.25 


6.6 




189 


254 


0.56 


0.70 


-0.17 


0.25 


0.25 


6.5 


quarter 


137 


137 


0.45 


0.46 


-0.14 





0.5 


6.4 


bridge 








0.30 


0.30 


-0.11 


0.25 


0.5 


6.3 




21 


21 


0.22 


0.22 


-0.09 


0.5 


0.5 


6.2 


hollow 


34 


42 


0.14 


0.14 


-0.07 
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Table 2 

Summary of results for the Br coverage dependence study. 



e 




^/N + B [meV] 


s [eA] 






[A] 


from DFT 


from DFT 


from DFT 


1/9 


8.91 


-2568 


0.31 


-0.12 


2/9 


6.3 


-2557 


0.32 


-0.12 


1/4 


5.94 


-2532 


0.30 


-0.12 


1/2 


4.2 


-2522 


0.27 


-0.13 


1 


2.97 


-1465 






2 


2.1 


5590 
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Fig. 1. A typical 3x3 supercell containing three metal layers, (light gray spheres) 
and one Br atom (dark gray sphere). The cell is 3(2.97) x 3(2.97) x 12.6 A in size. 
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Fig. 2. Electronic charge density, pd^), for the lowest energy configurations for 
several binding sites, summarized in Table 1. Few meaningful differences are seen 
outside of the adlayer region, even in a magnified view of the surface region (inset), 
e = 1/9. 
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Fig. 3. The charge transfer function, /S.p{z), for the lowest energy configurations 
for several binding sites, summarized in Table 1. The straight horizontal line is 

used to indicate Ap{z) = 0, and the vertical line indicates z = 5.4 A. The charge 
polarization in the surface region is easily observed. = 1/9. 
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Fig. 4. Four adlayer supercells. Small gray circles represent the top layer of Au(lOO) 
atoms, and the larger white circles represent adsorbed Br. Each cell has been re- 
peated in each lateral direction, and the repeated cells are shown by lighter shading, 
(a) 2 X 2 cell with 9 = 1/4, (b) 2 x 2 cell with 9 = 1/2, (c) 3 x 3 cell with 9 = 1/9, 
(d) 3 X 3 cell with Q = 2/9. The arrows indicate the adlayer lattice spacing, qbi, 
for each cell. 
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Fig. 5. Total energy from DFT calculations vs ZBr for several lateral positions. See 
Table 1 for more information. = 1/9. 
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Fig. 6. Sinusoidal corrugation potential, U {x, y). The grid lines indicate the bound- 
aries between surface unit cells, and the underlying Au(lOO) atoms are labeled. 
The grayscale image in the central square shows U{x,y), where black indicates the 
lowest value, meV, and white indicates U{x, y) = 137 meV. The gray shading is 
truncated above U{x,y) = 137 meV, and the contour lines indicate U{x,y) = 15, 
35, 55, and 137 meV as labeled in the plot. 
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2 4 6 8 10 

Fig. 7. Lateral interactions determined by DFT calculations along with an approx- 
imate analytic form. 
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+ 







8.91 A 



Fig. 8. A two-dimensional slice through the three-dimensional charge transfer func- 
tion, sliced through the two top sites to show the polar surface bonds associated 
with bridge site binding. The charge transfer function is for a 3 x 3 cell containing 
a single Br adsorbed at the bridge site. The circles indicate the adsorbed Br atom 
and the two Au atoms nearest to the Br, as labeled in the figure. The color scale 
indicates the local electron density, as indicated in the legend to the right. = 1/9. 
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1.3 




1 . 1 
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n 7 




















0.1 




-0.1 




-0.3 




-0.5 



O(x,y)=(0,0) a 
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0(x,yMl/4,l/4) a 
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Fig. 9. The surface dipole moment, s, vs. zbi for all considered lateral configurations 
with 9 = 1/9. It is quite clear, that even though the charge transfer function has a 
strong dependence on the lateral location, the surface dipole moment is essentially 
independent of lateral location and has a nearly linear dependence on zbt- The 
horizontal line indicates s = 0, and the heavy line indicates a linear fit to all data 
points. 
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Fig. 10. The surface dipole moment for four coverages, = 1/9 (a), 2/9 (b), 1/4 
(c), and 1/2 (d). The shaded symbols represent the configuration corresponding to 
the minimum energy for each coverage. 
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0.1 
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6 6.2 6.4 6.6 6.8 7 

ZBr[A] 

Fig. 13. The residence charge vs. Zbr for all considered lateral configurations with 
6 = 1/9. The linear dependence indicates that a charge/image charge description 
of the adlayer is not correct. 
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(I [meV] 

Fig. 14. The results of the off-lattice Monte Carlo simulations for an L = 32 
surface. Four different phases of the Br adlayer are seen in the simulations. The 
three ordered phases are seen as coverage plateaus. The insets show the short-range 
two point correlation function for each of these phases, shown in grayscale, where 
white represents places of highly coordinated atomic positions. An atom always sits 
at the origin of each of the insets, and "a" indicates the surface Au(lOO) lattice 
spacing, a. 
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Experiment, forward scan 

Experiment, reverse scan 

Simulation, L=32 




E [mV vs SCE] 

Fig. 15. The results of the off-lattice Monte Carlo simulations compared to the 
experimental cyclic voltammetry results. The experimental results are taken from 
Ref. [15] with permission from the authors. Note that peak heights are very difficult 
to calculate by Monte Carlo methods, and the peak positions are far more important. 
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